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Limiting behaviors of the Brownian motions on 

hyperbolic spaces 00 



H. MATSUMOTO (Nagoya) 



Abstract 

Using the explicit representations of the Brownian motions on the hy- 
perbolic spaces, we show that their almost sure convergence and the central 
limit theorems for the radial components as time tends to infinity are eas- 
ily obtained. We also give a straightforward strategy to obtain the explicit 
expressions for the limit distributions or the Poisson kernels. 



1. Introduction 

Hyperbolic spaces are non-compact Riemannian symmetric spaces of rank one. 
By classification, we have four types of hyperbolic spaces: the real one = 
SOo(l,n)/ SO(n), the complex one H™ = SU(l,n)/ SU(n), the quaternionic one 
Sp(l, n)/(Sp(l) x Sp(n)) and the Cayley hyperbolic plane. In this article we 
consider the limiting behaviors of the Brownian motions, that is, the diffusion 
processes generated by the Laplace-Beltrami operators on the first three types of 
the hyperbolic spaces. 

The hyperbolic spaces have negative bounded curvatures. The Brownian mo- 
tions on negatively curved manifolds have been studied in the connection of so- 
called the Liouville property by many authors and it is well known that the 
Brownian motions tends to infinity almost surely as time tends to infinity. See, 
e.g., Kifer |16j . Needless to say, the limit distributions are given by the Poisson 
kernels. 

On the other hand, the Brownian motions on the Riemannian symmetric 
spaces of non-compact type have been also studied by several authors since the 
work by Malliavin-Malliavin [20] ■ Among them, we refer to Babillot [3j, where 
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a central limit theorem for the radial components of the Brownian motions has 
been shown. 

The purpose of this article is to show that these properties are easily and 
directly shown for the Brownian motions on the hyperbolic spaces if we adopt 
the upper half space realizations of the hyperbolic spaces instead of the ball 
models. We can describe the same stories on the three types of spaces. At 
first, by solving the corresponding stochastic differential equations, we represent 
the Brownian motions in closed forms as Wiener functionals. Then, the almost 
sure convergence of them is readily seen from the representations. Moreover, by 
inserting the representations into the formulae for the distance functions, we can 
also show the central limit theorems for the radial components. 

For the computations of the limiting distributions or the Poisson kernels, we 
need some results on the distributions of the random variables defined by the 
perpetual (infinite) integrals in time of the usual geometric Brownian motions 
with negative drifts. The auxiliary results are given in the appendix and, by 
using them, we compute the Fourier transforms of the limiting distributions and 
the inverse transforms in direct ways. 



2. Real hyperbolic spaces 

For n ^ 1, let H™ +1 be the upper half space in R n+1 , 

H™ +1 = {z = {x,y) = (x 1 ,...,x n ,y);x e R n , y > 0}, 

endowed with the Riemannian metric ds 2 = y~ 2 (dx 2 + dy 2 ). The volume element 
is given by y~ n ~ 1 dxdy, the distance function d(z, z') is given by 

(2.i) cosh^z, z')) = \ x - x '\ 2 + y* + (y')\ 

in an obvious notation, where \x\ is the Euclidean norm. The Laplace-Beltrami 
operator is written as 

n d 2 d 2 d 

A r = y 2 \^ h y 2 (n — l)y — . 

^ dx 2 dy 2 dy' 

3=1 J 

For details of the fundamental objects on H™ +1 , see, e.g., Davies [6]. 

We first show an explicit expression as Wiener functional of the Brownian mo- 
tion on H™ +1 by soloving the corresponding stochastic differential equation. Let 
{W^ n+l \ B^ n+1 \ p( n+1 )) be the (n + l)-dimensional standard Wiener space with 
the canonical filtration {£>i n+1 ^} t > . Corresponding to the rectangular coordinate, 
we denote an element of W^ n+1 ' by 

(w(-), B(-)) or ( Wl (-) 7 ...,w n (-),B(.)), 

which is an R n+1 -valued continuous function on [0, oo) with Wi(0) = B(0) = 0. 
Then the Brownian motion on H™ +1 , the diffusion process with infinitesimal gen- 
erator A r /2, is obtained by solving the following stochastic differential equation 
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defined on (W^ n+1 \ B^ n+1 \ p(" +1 )) (see [32]): 

dXi(t)=Y(t)dwi(t), i = l,...,n, 

71 — 1 

dY(t) = Y(t)dB(t) —Y(t)dt. 

The unique solution Z z = z),Y(t, z))} t > , z — i x iV)i satisfying X(0) = x 

and Y(0) = y is given by 

Xi(t,z)=Xi + [ yexpiB^dwiis), 
Jo 

Y(t,z)=yexp(B { -^), 

where = B(s) - /is and p = n/2. {Y(t,z)} is a usual geometric Brown- 

ian motion with negative drift and it is easy to see that Z z (t) converges to the 
boundary as t — > oo almost surely. 

Now we consider the exponential functional A^~^ given by 

A { ~^ = [ exp(2B^)ds. 
Jo 

Then it is easy to see the identity in law 

(X(t,z),Y(t,z)) ( = w) (x + yw(4-%ye*p(Bt' t) )) 

for fixed t > 0. 

An explicit expression for the density of the distribution of (A\~^\ b[~^) is 
known by Yor [27J and, by using it, Gruet [12] has shown an expression for the 
heat kernel of the semigroup generated by A P . For the classical expression, see 
Davies [6]. We also refer to [2j [H EEJ E3] for related topics. 

We combine the identity in law with formula (12. IB . Then we get 

«*h(d(Z(M),*)) ( = W) i{|^(A ( " M) )r + l}ex P (-^) + iex P (^). 

Since A[ ^ converges as t —>■ oo and log(cosh(u)) = u • (1 + o(l)) as u — > oo, we 
readily get the following central limit theorem. 

Theorem 2.1. The probability distribution of 'y/t 1 (d(Z(t, z), z)—nt/2) converges 
weakly as t — > oo to the standard normal distribution. 

Recall formula A r d(z , ■) = n coth<i(zo, •)■ Then, by Ito's formula, we get 



n „ 

d(Z z (t),z)=J2 / 
Jo 



1 X 'M " * iw i 



^J sinh d(z,Z z (s)) y 



smhd(z,Z z (s))\ y 

Th f 

+ — / coth d(z, Z z (s))ds, 
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from which the theorem may also be proven. 

Next we recall Dufresne's identity (Theorem IA.1I in the Appendix) in law 
^(-/ t ) (law ^^J- 1 f or a Gamma random variable 7^ with parameter /1. Then, for 
a bounded continuous function ip on R", we obtain 



E[ V (X(t,z))}=E[ V (x + yw{A { r /2) ))] 



l_^M c -«dt I <p{x + V ) fn _J l ^ e X p[--S 7K[ )dr) 
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o r(n/2) 7 R ^ V " (27q/72^/2 *\ % 2 /2t 

^ Vo r(n/2)WV V 2/ 2 / 

¥>(OPn+l(f -Z,2/)df, 



R' 



R» 

where 

_ 2- 1 r((n + l)/2) y" 

7r(n+l)/2 ( y 2 + |£|2)n' « ^ ' 

and we have used the duplication formula for the Gamma function. 
Hence we have proved the following. 

Theorem 2.2. For any (x,y) G H" +1 , X(t,z) converges almost surely as t — > 00 
and the density of the limit distribution is given by the Poisson kernel p n+ i(£ — 
x,y). In particular, when n = 1, the limit distribution is Cauchy. 

We end this section with mentioning on the Poisson kernel in the Euclidean 
spaces and on the Fourier transforms. The Poisson kernel on R n+1 of the hyper- 
plane {y = 0} is given by 

r((n + l)/2) y 



n (n+l)/2 ^2 _|_ |£|2)(n+l)/2' 

which is different from p n +i{£,,y) for n = 2. The Brownian motion on the hyper- 
bolic plane H 2 is a time change of the 2-dimensional standard Brownian motion, 
and the Poisson kernels coincide. 

It is well known that the Fourier transform of q n +i(£, y) in £ is the simple 
exponential function, 

For the hyperbolic spaces, we can show, for example, 

^{\;y) = [ e^ T ^p 3 (^y)d^ = y\X\K 1 (y\X\) 
7r2 

and 

M*;y)= [ e^ T ^ PA (^y)d^=(y\x\ + i) e -y\ x \, 

JR3 
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where K\ is the modified Bessel function. By virtue of the strong Markov prop- 
erty, we can easily show that the distribution of X{r a ) for the first hitting time 
r a of the Brownian motion {Z z (t,z)} at the level y = a, a > is determined by 
the characteristic function written by 

£[exp(v^T(A,X(r a )»] = e^^J, A G R" 

(p n (A;a) 

3. Complex hyperbolic spaces 

Let H™, n ^ 2, be the upper half space of C n given by 

{z = (z h z 2 , z n ) = (z 1 ,z)e C n ; h(z) = Im(>i) - \z~\ 2 > 0}, 
endowed with the Bergmann metric 

n 

ds 2 = - d z .dg^log(h)dzjdzk. 
j,k=i 

The unit ball {\z\ < 1} in C n with the Bergmann metric 

n 

- <9^<3^1og(l - |2| 2 )d^d^. 

is isometric with H™. For details, we refer to pO, [8], QUI El]- We should be aware 
of difference of conventions. The curvatures of these manifolds are bounded and 
negative, but are not constant (cf. p. 190, [Tu]). 

We change the first coordinate by X\ = Re(zi)/2 and y = h(z) 1 ^ 2 . Then we 
have the same realization of the complex hyperbolic space SU(l,n)/ SU(n) as in 
Venkov [24|: if we write Zk = x k + y—lyu, k = 2, n, the Riemannian metric is 
written as 

1 1 n 1 / - \ 2 

ds 2 = —dy 2 + — ^2(dx 2 k + dy 2 k ) + — ( dxi + ^{x k dy k - y k dx k ) J , 

and the distance function d(z, z') is given by 

(cosh(d(^))) = i -^- 7 p , 

where 

n 

(3.1) ® = y 2 + \z'-z\ 2 and ip = x[ - x 1 + ^{y' k x k - x' k y k ). 

k=2 

The Laplace-Beltrami operator is written by 

d 2 d 2 d 
dx\ ay 1 ay 

(3.2) „ , , d a \2 / a a 



k=2 k V 



n +Vk^~ + ~ x k 

dx k axx J \dy k dxi 
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Letting (W^ 2n \ B^ 2n \ p( 2n )) be the (2n)-dimensional standard Wiener space 
and denoting an element of W^ 2 ^ by 

(B(-),w 2 (-),w 3 (-),...,w 2n (-)) or (B(-),w 2 (-), $(■)), 

we can check that the Brownian motion {Z(t)} on H™ with Z(0) = (x±,y,z 2 , 
z n ), z k being identified with (x k ,y k ), is given by 



X(t) = xx 



+ / Y(s) 2 dw 2 (s) + 2j2S k (t) 

J ° k=2 



(3.3) 



Y(t) = yexp(B(t) - nt) 

X k (t) =x k + Y(s)dw 2 k-i(s), 
Jo 

Y k (t) = y k + Y(s)dw 2k (s), k = 2,...,n 



where we have used the trivial notations X(t),Y(t), Xk(t),Yk(t) for the compo- 
nents of Z(t) and S k (t) is the stochastic area enclosed by {(X k (s), Y k (s))} < s < t 
and its chord, 

S k (t) = l [ (Y k (s)dX k (s) - X k (s)dY k (s)). 
* Jo 

{Y(t)} is again a usual geometric Brownian motion with negative drift and 
Z(t) converges as t — > oo. Hence we easily obtain the following central limit 
theorem. 



Theorem 3.1. For the Brownian motion {Z(t)} on the n- dimensional complex 
hyperbolic space, the probability law of \fi (<i(Z(t), Z(0)) — nt) converges weakly 
as t — ► oo to the standard normal distribution. 

Next we compute the limiting distribution of (X(t), Z(t)) as t — > oo, where 
Z(t) = (X 2 (t),Y 2 (t), ...,X n (t),Y n (t)). 

If we consider the ball model, we obtain the Poisson kernels as the densities of the 
image measures of the uniform measure on the sphere by the isometries. However, 
since the same strategy works in the more complicated case of the quaternionic 
hyperbolic space whose geometry has not been well understood (see a recent 
work by Kim-Parker [TTj and references cited therein), we give the following 
straightforward computations. 

For this purpose we first fix t and consider the characteristic function. As in 
the previous section, we set 
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B [ s M) = B(s)-fis and \i = n. For the stochastic analysis on H™ and H™, we need 
to consider these two exponential functionals. Then, by the expression (I3.3p . it 
is easy to see that, for fixed t > 0, {X(t), Z(t)) is identical in law with 

where ^2 denotes the sum over k = 2, n, <p(t) = J2(yk w 2k-i(t) — %kW2k{t)) an d 



Sk{t) = - I (w 2 k(s)dw 2 k~i(s) - w 2 k-i(s)dw 2 k{s)). 



Hence we may write, for any bounded continuous function g on R 2 ( n ^, 



E[e^ x Vg(Z(t))} = E 



x E 



L fc=2 



{B(s)},w{A^) 



Then, applying the Levy formula for the characteristic function of the stochastic 
area (cf. [15], p. 473), we get 



E[ e J =T * x ®g(z + Z(t))] = E 



py A\ w x 



sinh(p?/ 2 y4j 



2 A {-ti Ai-^h\\E A t^)l 



exp (l-py 2 A\ w coth(py 2 Aj 



2A 



(-m) 



Moreover we take the conditional expectation given {-B(s)} s > or {^(s)} s > to 
obtain 



E[e^ pX Vg(Z(t))] = e^ px E 



-P 2 V 4 A^/2 



P 



R2(n-1) 
n-1 



27rsinh(p|/ 2 y4 



HO* 



-pcoth(p y 2 A(-^)|C| 2 /2 d ^ d7? 



where C = (£, ^) = (6, £n, ^n)- 

Now we put, for q = (q 2 , q n ), r = (r 2 , r n ) G R n_1 , 

( 7 (C) = exp(v /Z T((q,0 + (r^))). 
Then, earring out the Gaussian integral with respect to £ and 77, we get 



E 



expi >/=T(pX(t) + ^(g fc X fc (t) + r k Y k {t))) 



-if E 



cosh(py 2 A[ ^) 



n-l 



-Ftanh(py 2 A< m) ) 



s 
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where 



/ = fip, q, r) = px + ^2(q k x k 



TkVk) 



and 



7-1 7-1/ v (g* + M/fc) 2 + (rjfc - par*) 
F = F(p, q, r) = 2^ — 



Now, letting t — > 00, we obtain the following. 



Proposition 3.2. For any p G R, q, r G R n , one /ias 



(3.4) lim F 



expi >/=T(pX(t) + ^(g fc X fc (i) + r k Y k (t))) 



if E 



-p 2 y 4 A ( ~ n) /2 



cosh(py 2 A 



(-nh 



n-l 



-Ft a nh(py 2 A { - n) ) 



Denote the right hand side of H3.4[) by I(p, q, r). By using the joint Laplace 
transform of A^ n ^ and A^, given by Corollary IA.5I in the appendix, we obtain 



Hp, q> r ) 



-1/ 



1/ 



-pV^ n) /2 



At n) 



u 



n-l 



cosh(pn 2 n) 

X e -F^HPV 2 u)p^ A (-n) e 



py 



o 2 n r(n) \sinh(pn 2 n) 



n+l 



-py 2 coth(j«/ 2 u)/2 



1 



n-l 



cosh(p?/ 2 n 

Then, changing the variable, we see that I(p, q, r) is equal to 

n+l 



e -Ftanh(py 2 u)^ 



U. 



f-i f {py 2 Y ' 



1 



2 n Y{n) J \ v sinh(n) J \cosh(n 



■ \ n— 1 

» e -py 2 coth(u)/2 



if p > and to 



f-Lf {-py 



2\n r°° 



x exp , - £ (9fc+ ^ )2 + (rfc - pXfc)2 tanh(n) ) dn 



, \ n+l / 1 \ n— 1 
1 ' ' 1 ' e py 2 coth(u)/2 



2 n T(n) J \sinh(n) / \cosh(w 



x exp , £ (9fc+ ^ )2 + (rfc - pXfc)2 tanh(n)) dn 
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if p < 0. From these expressions, we can take the Fourier inversion 

f n (x', 7; z) = 1 1 / I(p, q, r)e-^ x ' + ^< +r ^dpdq 2 ■ ■ ■ dr r , 
For the integral with respect to q k when p > 0, we note as usual 

(qk + pyk) 2 



— tanh(u) + V-Lq k (x k - x' k ) 

2p 

tanh(it) , , — - . . , ,x2 

■{qk+pyk - v-ip{x k - x k ) coth(u)j 



2p 



^lpy k {x k - x' k ) ~ 2 ( x k - x 'kf coth(w). 



We do the same computations also for the other variables and for p < 0. Then, 
after some manipulations, we obtain 

f n {x , Z , z) 

„2n f poo / i \ 2n 

/ {p^e^^dp / (-^-- e -^ coth ^ 2 du 



p 2 - 1 cos(2 W )dW ( — * ) e -««>du, 



where we have made a simple change of variable for the second equality and <p 
and $ are given by (13. ip . 

For the second integral, we change the variable by k = coth(w) to obtain 

o„ 2n poo poo 

f n (x', 7; z) = -J— / p 2 "- 1 cos(2^)dp / e~^ k (k 2 - ly^dk. 
71-1 W J J 1 

Now we recall the following integral representation of the modified Bessel function 
(cf. Lebedev [H] p.119 or [llj p.322) 

Then we obtain 

on+1/2 2n poo 

fnW,7;z) = r + 1/t £_ 1/3 J Q P n - 1/2 ■cos(2<pp)K n - 1/2 (*p)dp 
For the integral on the right hand side, we may apply the formulae 

I™ x x K,(ax) cos(fe)dx = 2 x - 1 a - x - 1 T{ fX + * + - )r( - + * ~ M ) 
Jo 2 2 



, /x + A+l 1 + A-/X 1 _b_ 
[ 2 ' 2 '2' a' 
(cf. [llj p. 747) and F(n, a, a; z) — (1 — z) _n . Then we obtain 

(3.6) / n (x , 2 ; z) = — — iF (n; -— ) 



7r n<J,2n 1 U V ' <|)2 / 7T n (4<£ 2 + <3> 2 ) n ' 



10 



H.MATSUMOTO 



Theorem 3.3 (cf. [7j). For any z G H™, (X(t),Z(t)) converges almost surely as 
t — > oo and the density of the limit distribution on R 2n_1 is the Poisson kernel 
given by (13.61) . 

4. Quaternionic hyperbolic spaces 

For the quaternion hyperbolic space Sp(l,n)/(Sp(l) x Sp(n)), n ^ 2, we follow 
the conventions in Venkov [21]. See also Helgason [13J, Lohoue-Rychner [19J, and 
Kim-Parker [T7] for the basic properties. For n ^ 2, let H" be the upper half 
space in C 2n , 

H™ = {z = (z 1 ,z 2 ,...,z 2n ) = ( Zl ,z) G C 2 ";Im(^) > 0}, 
with the Riemannian metric 

dy 2 1 n 
ds 2 = — + — V7dz fc dz^ + dz n+k d~z~^) 

y y t^2 

i / n \ 2 

+ — \^dx 1 + Im Y^(z^dz k + z—^dz n+k ) J 

n 

dz n +i + ^^(^n+fcdz fc — z k dz. 



k=2 

1 

+ 7 



fc=2 



where Z\ — X\ + y— Ty. We will write = x k + y/ — iy k for k = 2, 2n. Note 
that the first and (n + l)-th components, Zi and z n+ i play special roles. 

The volume element is y~ 4n ~ 3 dxidy n&=2 ^ x k^Vk an d the distance function 
d(z, z') is given by 

,4,, {a ^ m ^ m±^^SL , 

where 

n 
fc=2 

n 

<Pi = x i -xi + ^2((y' k x k - x' k y k ) + (y; +fc x n+fc - x' n+k y n+k )), 
(4-2) fc=2 n 

<^2 = x' n+1 - x„+i + ^2((x' k x n+k - x'„ +fc a; fc ) + - y' k y n+k )), 

k=2 

n 

<^ 3 = y'n+i ~ yn+i + ^2((x' k y n +k - y' n +k x k) + (y'k x n+k - x' n+k y k )). 

k=2 

Note that (p^s do not depend on y. 
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The Laplace-Beltrami operator A q may be written in a convenient way as 
d 2 d 2 d ( d 2 d 2 

A '=^ + ^- (4 " +1) ^ +!/1 fe + ^ 



k=2 



( 4 - 3 ) + h^-x k ^ + y n+k 



d k d d d 

dx k dx x ox n+ i dy n+ i 

d u d d d x " 



dy fe dxi dx n+1 dy n+l 

d d d d 



/ d 9 d d \ 2 

\dy n+k dxi dx n+ i oy n+1 J 

Note that the coefficients of d 2 / dxidx n+ i, d 2 / dxidy n+ i, d 2 / dx n+ idy n+ i are zero. 
We can describe the same story as for the complex hyperbolic space if we consider 
a 4 x 4 skew-symmetric matrix instead of 2-simensional one. 

At first we give an explcit expression for the Brownian motion, the diffu- 
sion process with generator A 9 /2, on H™. Let (W^ n \ £>( 4n ), p( 4n )) be the tri- 
dimensional Wiener space and denote an element in W^ 4ra ) by 

(B 1 (-),B(-),W 2 ,l(-),W 2 ,2(-), -,Wn,l(-),Wn,2(-), 

B 2 (-), B 3 (-),w n+2:1 (-),w n+2j2 (-), w 2n ,i(-), w 2n , 2 (-)). 

Then we can check that the Brownian motion (X(t),Y(t), Z(t)) starting from 
(x 1: y,z) is given by 



t 

2, 



X 1 (t)=x 1 + / Y(sYdB 1 (s) 
Jo 

n nt 

+ J2 / {Yk(s)dX k (s) - X k (s)dY k (s) 

k=2 Jo 



+ Y n+k (s)dX n+k (s) - X n+k (s)dY n+k (s)}, 
Y(t)=yexp(B(t)-(2n + l)t), 

X k (t) = x k + / Y(s)dw kA (s), 



Yk(t)=yk+ Y(s)dw k>2 (s), k = 2,...,n, 

Jo 

X n+1 (t) =x n+1 + [ Y(s) 2 dB 2 (s) 
Jo 

n ft 

+ J2 {-X n+k (s)dX k (s) + X k (s)dX n+k (s) 

k=2 J ° 

+ Y n+k (s)dY k (s) - Y k (s)dY n+k (s)}, 

Y n+1 (t)=y n+1 + ! Y(s) 2 dB 3 (s) 
Jo 
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n pt 

+ Y1 / {-Yn +k (s)dX k (s) + X k (s)dY n+k (s) 

k=2 J ° 

- X n+k (s)dY k (s) + Y k (s)dX n+k (s)} 

X n+k (t) = x n+k + / Y(s)dw n+ktl (s), 
Jo 

Y n +k{t) = y n +k + Y(s)dw n+k , 2 {s), k = 2,...,n. 

Jo 

Then, from (14. ip . it is easy to show the following central limit theorem. 



Theorem 4.1. The probability law of (d(Z(t), Z(0)) — (2n + \)t)/yi converges 
weakly as t — » oo to the standard normal distribution. 

Next we show that (X(t),Z(t)) converges in law as t —>■ oo. To identify the 
limit distribution, we set 



f n (x,z,z) 



)4n+l 



T(2n) 



y 



2(2n+l) 



7T 



2n 



($ 2 + 4(^ + ^I + ^)) 2 " +1 ' 



where $ and (p^s are given by ( 14. 2p . / n is the Poisson kernel of the boundary 
dH£ = {y = 0}. 

Theorem 4.2. Z(t)), valued in R x C 2 ^ -1 -*, converges almost surely as 

t — > oo and £/ie density of the limit distribution is given by f n (x', z 1 ; z). 



In the following we give a proof of Theorem 14.21 At first we consider the 
characteristic function of (Xi(t), Z(t)) for fixed t. For convenience we put 



X° k (t) = X k (t) - x k 
and 



Y(s)dw k , 1 (s), Y k °(t) = Y k (t)-y k = / Y(s)dw k , 2 (s) 
o Jo 









( x k (t) \ 




( x° k (t) \ 


Ok — 


Vk 

Xn+k 

\y n +k ) 




Y k (t) 
X n+k {t) 

\Yn + k(t)J 




Y k ° 

^~n+k 

\Y° +k (t)/ 



Moreover, £ = *(^i,6,6) e R 3 ,w k = t {u k ,v k ,u n+k ,v Tl+k ) G R 4 , w = (w 2 , w n ), 
we set 

(t) = + 6^n+l(0 + (*), 

c/ fc (t) = (iw fc , e fe (t)>, c/ fe °(t) = k, e°(t)>. 

We throughout denote by t Q the transpose of a matrix Q. Then the characteristic 
function is 



<p(t) = E[exp{V=l(C 1 X 1 (t) + ^X n+1 (t) + 
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+ \/-L^2(u k X k (t) + v k Y k {t) + u n+k X n+k {t) + v n+k Y n+k {t)) }] 

= £[exp(v^(*(0 + 

where the summation is taken over k = 2, ...,n. 

To compute the characteristic function, we introduce a 4 x 4 skew symmetric 
matrix S given by 



( 


6 


-6 


-&\ 











6 











-6 


-6 


o / 



Then we have 

*(f) + 5>*(*) 



^ + / F(s) 2 (eid J B 1 ( S ) + &dfl 2 (s) + e 3 dS 3 (s)) 
./o 

+ £<s0 fe + Wjfc , e2(f)> + £ Ase2( a ), de°( s )), 

Jo 



where ip = £ 1 x 1 + £ 2 x n +i + fsj/n+i + E( w fcA)- Note that {Ej=i £j 5 j( s )} is 
identical in law with |f | = fe 2 + £ 2 2 + £ 3 2 ) 1/2 - 

The eigenvalues of S are ±\/— 1|£| and the multiplicities are two. Moreover 
there exists an orthogonal matrix Q such that t Q'BQ = K is of the standard form. 
We take 



/ 



Q 



& 
u id 

1 

u m 



i*i 
o 








\ 



and K 



/o o \ 

ICI o 

o o o -lei 

\o o lei o / 



We also put = t Qw k) 6 k = l Q9 k and 

©2(0 = (^(t)S k \t)J° n+k (t)XlM) = 'q§2(0. 

By the rotation invariance of the probability law of Brownian motions, we see that 
{@2(s)} is a simple time change of a 4-dimensional standard Brownian motion. 
Under these notations, we have 



(Z6 k + w k , 92(f)) = (K9 k + w k , 92(f)) 



and 



(4.4) 



f(~92( S ),d92( S )> 
Jo 



= lei /V*°(«)d^(«) - Y k °(s)dXl(s) 
Jo 

+ ^ +fc («)di; +fc ( S ) - t° +fe («)d^ +fe ( S )}. 
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As in the previous sections, we set 



A 



(-*0 



exp(2B { -^)ds and A\ 



E>i ^ = B(s) — us and /i — 2n + 1. Then, by taking the conditional expectation 
given and 6°(t), k = 2, ...,n, and applying the Levy formula, we obtain 



<p(t) = e^^E 



x 



n^[ex P fv=T Ase2( a ),de2( a )> 



-l^(^ fc + «; fc ,yV? fc (4^ ) )) 



2j(-c) \2(n-l) 



smh(|e|2/ 2 ^-^^ 



ex P ( (l-ieb^cothdeb 2 ^)) 



2 



(-^)M2 



2„,4 T(-A*) 



ex P --ieiv^ 



k c l 



R4 (n-i) \27rsinh(|£|?/ 2 v4 



2(n-l) 



exp/V^l^fW + ™*,G> - ~|£| coth(|e|y% ( - M) )|C fe | 2 )") dC 2 • • - dCn 



where {(PV^s), W„(s))} is a 4(n — l)-dimensional standard Brownian motion, 
independent of {y(s)} or {B(s)}. We carry out the Gaussian integral over R 4(n_1 ) 
to obtain 



<p{t) 



-it E 



cosh(\£\y 2 A 



ho* 



2(n-l) 



-|e| 2 2/ 4 4" /l) /2-i ? tanh(|e|j/ 2 4-' i) )/2|e| 



where F — \K9k + Wk\ 2 - Hence, letting t tend to oo, we obtain 



lim ip(t) 

t—>oo 



-11, E 



cosh(|£|?/ 2 A 



Hh 



2(n-l) 



-|e| 2 2/ 4 ^t /l) /2-i ;, tanh(|e|j/ 2 yl^ w )/2|e| 



2 



Now, applying ( 1A.9I) . we obtain the explicit expression for the Fourier trans- 
form of the Poisson kernel f n . 



Proposition 4.3. Under the notations above, the Fourier transform of the limit 
distribution of (X (t) , Z (t)) as t —>■ oo, is given by 
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e ,/=lil> /\£\y2\2n+l roc / j \ 2(n-l) / 1 \ 2(n+l) 



2 2n+1 r(2n+l) 7 \cosh(u)/ Vsinh(M 



x exp^-^|£|?/ 2 coth(w) - ^|tanh(w)^dw. 
We invert the Fourier transform. That is, setting x' = (x , i,x' n+1 ,y , n+1 ), 6' k - 

t ( X 'kiy'ki X 'n+k^y'n+k) anC ^ 

we compute 
f n {x , Z j z^j 

i 4n— 1 

f*oo fU\y2\2n+l / i \ 2(n-l) / 1 \ 2(n+l) 

X 



2 2n+1 T(2n + 1) \cosh(u) / \sinh(M) 
x exp^-i|e|y 2 coth(n) - t ^^± ^ \K6 k w k \ 2 j du. 



Recall the definitions, = *Qwfc and = 'Q^fe- Then, changing the order of 
the integrations, we have 

f n (x i z j z) 

2(2n+l) /-oo / \ 2(n-l) / -, n 2(n+l) 



(27r) 4 "- 1 2 2n+1 r(2n + 1) 7 V cosh 0)/ V. sin M M 

x /" e ^($,a:-a ; ') e -|$|2/ 2 coth(w)/2|^|2 n +l d ^ 

n 

^E^fc ,? k -^>-tanh(«) E \K9 k +w k \ym JJ ^ 

fc=2 



R4(n-1) 



We can easily carry out the third Gaussian integral since Q G 0(4) and we obtain 

, ~, y2(2n+l) roc / j_ \ 4n 

f«tf^'^ = (47r)2«+ir(2n + l) 7 VsTnhHj ^ 

X / e V^{<€.»-a/>+|C[0(S;fl > )}-[€[*coth(u)/2|^|4i»-l d ^ j 

where $ = y 2 + £ |0 fc - = y 2 + £ |0 fc - 0£| 2 and 

0(#> = (-K"^fc, A) = ^{y'kXk - x' k Vk + y'n+k^n+k - x'n+kVn+k)- 

For the right hand side, changing the variables by k = coth(w), we first com- 
pute the integral in u. Then, by using formula (13.51) again, we get 

i \ 4n poo 
——-) e -m^ H u)/2 du= / -2 _ 1) 2n-l e -[ € |**/2 dA . 

sinh(tt) J J x 

Y{2n)( 4 \ 2n - l! \ (m 
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Moreover, by definitions, we see 

(t,a/-x) + \Z\<l>{e,e > ) = -(t,<p), 

where (p — t (cp l , <p 2 , ^3) is given by 04.21) . 
Combining these identities, we obtain 



f n (x , z , z) 

*-T7J L e ~^" )K ^ (flfl) l«| 2 "" V2d « 



y2(2n+l) 



8(2n + l)7r 2n + 3 / 2 $ 2 

^/ R3 ^^ l6 ^-i/ 2 (fiei)iei 2 "- 1/2 de- 



y2(2n+l) 



8(2n + l) 7r 2n+3/2 < | ) 2n-l 

Moreover, changing the variables by the spherical coordinate, we obtain 

4 vr y2(2n+l) 



f n {x,z,z) 



\(2n + 1)^+3/2^-1/2^1 

x J™ \ 2n+1/2 K 2n _ l/2 (^r^ sin(|^|r)dr. 



For the integral on the right hand side, the following formula is available (cf. [H 
p.747): 

I™ x x K^ax)sin{bx)<\x = 2V A - 2 6 x r ( 2 + ^ + A ) r ( 2 + A ~ ^ ) 
Jo 2 2 

+ p + A 2 + A-ji3 b 2 

b ( 

In our case A = 2n + 1/2, = 2n — 1/2 and 



v 2 ' 2 '2' a 2 ' 



( 2 + n + \ 2 + A-^ 3 / fc 2 ^ ~ (2n+1) 

^ 2 ' 2 '2' a 2 ^ 1 V a 2 
Hence we may apply this identity and we arrive at our result 

2 4n+1 T(2n) y 2 ( 2 "+!) 



fn Z , Z) 



vr 2n ($2 _|_ 4|^,|2)2n+r 



A. Perpetual integrals of geometric Brownian motion 

In this appendix we consider two perpetual integrals of geometric Brownian mo- 
tions. Let B = {B(t)} t > be a one- dimensional Brownian motion with Bq = 

defined on a probability space {tt,F,P). For fi > 0, we set = {B { t m) = 

B(t) — fit}, a Brownian motion with negative constant drift — fi. Then Dufresne's 
perpetual integral is defined by 



(A.l) 



A^-fi = / exp(2.B^)ds. 
Jo 
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We also consider another integral 

Jo 

Then the following is known: 

Theorem A.l (Dufresne [9]). Let'j^ be a gamma random variable whose density 
is given by(T(/j,))- 1 x' M - 1 e- x . Then is distribued as (27^) 1 and, accordingly, 



Remark A.l. Several different proofs of this theorem are known. In particular, see 
Yor [26J. The density of the exponential functional A t = J Q * exp(2B s )ds for fixed t 
has been obtained by Yor [27] and the joint distribution of (At, a t ) in an obvious 
notation has been studied in [TJ. See also (22J, [23J, [25] for several results and 
applications of these perpetual integrals and exponential functionals. Recently 
Baudoin-O'Connell [5] has shown several formulae, including (IA.20 below, for 
the exponential type Wiener and discussed their close relation to the theory of 
quantum Toda lattice. 



What we need in Sections 4 and 5 is the following explicit expression for the 
iditional Laplace 
by Yor [2E]. We set 



conditional Laplace transform of Aio^ given aL M \ which was originally obtained 



9 2 M 

which is the density of the random variable aL^ or 2/7^. 
Theorem A.2. For A > and v > 0, it holds that 



E 

(A.2) 



exp 



1 / a' n 2 "" j 



h(v) 



2T(2/j,) V.sinh(Av/2) 



exp —A coth 



(*))■ 



We have this nice result only for the particular choice of A t and at, that is, it 
is available only when the ratio of the coefficients in the exponential functionals 
is two. 

We give another proof of the theorem for completeness. Note that by letting 
A tend to in (|A.2j) . we obtain Theorem IA.11 

For this purpose, we consider the Brownian motion {B^ = B t + /it} with 
the opposite positive drift and set X x (s) = xexp(Bs), which defines a diffusion 
process with infinitesimal generator 

1 2 d 2 (\ \d 



2 s d^ + U + / Td*- 
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Letting r 2 be the first hitting time of {X x (s)} at z, we set for A > and kgR 



VAX 



E 



exp(-y ^ X x (s)- 2 ds + Xk ^ X x (s)- X te\ 



In [22] we have considered the case of k = and showed that v z (x) may be 
represented by means of the modified Bessel function to give another proof of 
Theorem IA.1I Following the same line, we first show a representation for v z (x) 
by means of the Whittaker function. 

Let W K ^ be a Whittaker function: if // — k + 1/2 > 0, 



(A.3) 



-2/2 M +l/2 poo 

r(//-K + i/2) Jo 



From this expression it is easy to see lim 2 ^ oc = when \k\ is small. We 

also recall that W Ksft solves the equation 



W"(z) 



k ^ - (1/4) 



4 z 



W{z) = 0. 



Proposition A.3. For pi > 0, A > and k G R, it holds that 

»- 1/2 W^(2\/x) 



(A.4) 



VAX) 



W K ^(2X/z)- 



Proof. We have only to consider the case of k < 0. The general case can be 
shown from this case by analytic continuation in k. Note that, if k < 0, v z (x) is 
monotone increasing in x(> z). 

At first we note that vAx) is a solution for 



x v"(x) + ( - + fi \ xv\x) 



( 'A 2 



- 



\2a; 2 



Ak 



x 



v{x) 



and satisfies 
(A.5) 



v g (x) 



and 



lim t>*(x) 



0. 



We now change the variable by £ = A/x and set 

Then, by straightforward computations, we see that (f) satisfies 

2k pi 2 -(1 /4)" 



4>(() = 0. 



By considering the boundary conditions ( 1A.5I) . we can easily show 



0(0 



W^(2£) 

a; w«„(2a/z) 



and hence the result ( 1A.4I) . 



□ 
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Proposition A. 4. For jj, > 0, A > and k 6 R, it holds that 



(A.6) E 



exp 



--\ 2 A^ + A/sa^O 



r(/x-« + 1/2) 
r(2//) 



(2A)^w ( 2A ). 



Proof. By the symmetry of the probability law of Brownian motion, {— B t } 
{B t }, we have 



(law) 



lim v z (l) = E 

z— >oo 

= E 



<>xp( --A 2 
<>xp( -^A 2 



-2B 



(m) 



ds + Ak 



,2B 



(-*») 



ds + Ak 



e~ Bs ds 



r(-m) , 

e Bs ds 



On the other hand, by using the fact on the Whittaker function 



(A.7) 



W K ,Jz) 



r(2/i) 



t(ji-k + i/2) 

we see from the expression ( 1A.4I) 

r(^-« + i/2) 



lim u z (l) 



r(2 M ) 



(2A)^ 1 / 2 ^ K , M (2A). 



□ 



Remark A. 2. The asymptotic behavior OA .70 of can be easily shown by the 
definition of the Whittaker functions. 

Now we are in a position to complete our proof of Theorem IA.21 By ( 1A.3I) 
and (1A.6I) . we have 



E 



exp 



-x 



r(2/i) 



(2A) 



2/j, 



-2At^/i-K-l/2/2 _|_ ^-j^+k-1/2^ 



Now we change the variable by e Xv = 1 + 1 1 or t = (e — 1) 1 . Then some 
elementary computations show that this integral is equal to 



X 2 M 



1 roc 



-Acoth(Au/2) 



2r(2/i) 

This completes our proof since 



sinh(Av/2) 



2|t+l 



e A ™d?;. 



exp 



E 



exp 



fi(v)e XKV dv. 



□ 
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Corollary A. 5. Define another perpetual integral A^o by 

POO 

(A.8) = / exp(4S(-"))d S 



o 



and let /2(f) be the density of Aio^ or (27^) 1 . Then one has 



E 

(A.9) 
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